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INTERACTION  BETWEEN  A  SHOCK  WAVE  IN  AN  ELASTO- PLASTIC 


MEDIUM  AND  A  NONRIQID  WALL 


G.  M.  Lyakhov  and  N.  I.  Polyakova,  Moscow 


The  interaction  between  a  wall  and  a  wave  in  an 
elasto-plastic  medium  with  an  alternating  sign  of 
curvature  of  the  stress-strain  dependence  a  =  a(e) 
was  examined  earlier  [ 1 ] .  The  case  of  small 
stresses  lying  on  the  concave  section  of  the  dia¬ 
gram  a  =  o(e),  when  the  wave  ha3  no  pressure  jump 
at  the  front,  has  been  investigated.  Below,  based 
on  earlier  studies  [1,  2],  we  will  give  a  solution 
to  the  problem  of  the  interaction  between  a  plane 
wave  and  a  nonrigid  wall  or  the  boundary  of  media 
in  the  case  of  large  stresses,  corresponding  to 
the  section  of  the  diagram  o(e)  which  is  convex 
relative  to  axis  e.  Here  the  wave  is  a  shock  wave. 


1.  We  will  examine  a  model  of  a  medium  with  nonlinear  dependence 

a  =»  a(e).  With  small  values  of  c  we  have  d2o/de2  <  0  and  with  the 

large  values,  d2o/de2  >  0.  Compression  and  unloading  when  a  <  a  occur 

s 

elastically,  i.e.,  according  to  one  law,  but  when  a  >  a_  they  occur  ac- 
cording  to  different  laws.  Secondary  loading  takes  place  according  to 
the  law  of  unloading  up  to  the  stress  achieved  on  first  compression. 
Such  a  model  is  applicable  to  the  ground  and  certain  solid  media. 

We  will  examine  large  stresses,  therefore  we  can  consider  the 
elastic  section  as  small  and  approximate  the  relation  a  =  o(e)  under 
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load  with  two  straight  lines  (Fig.  l). 


We  will  examine  two  variants  of  the  approximation  of  the  unload¬ 
ing  curve.  In  the  first  case  we  will  assume  that  unloading  takes 
place  at  a  constant  residual  deformation,  i.e.,  along  a  line  parallel 
to  axis  a,  in  the  second  case  the  unloading  line  is  parallel  to  the 
second  section  of  the  linearized  compression  diagram  a  »  a(e). 

In  Fig.  la  and  b,  broken  line  1  corresponds  to  the  loading  line, 
straight  lines  2  to  the  unloading  lines  under  constant  residual  defor¬ 
mation,  and  straight  line  5  to  the  unloading  lines  parallel  to  the 
second  section  of  broken  line  1. 

Wave  propagation  in  the  first  law  of  unloading  has  been  examined 
by  V.  N.  Rodionov,  A.  Ya.  Sagomonyan,  and  other  authors;  the  reflec¬ 
tion  from  a  rigid  wall  for  the  approximation  of  compression  curve 
a  =  a(e)  by  a  straight  line  was  examined  by  S.  Kallskl  and  Ya.  Osieski 

[?3. 

We  will  change  from  parameters  o  and  e  to  the  parameters,  pres¬ 
sure  and  volume  V. 

For  the  propagation  of  a  plane  wave  in  a  medium  unbounded  in  a 
direction  perpendicular  to  its  movement,  or  in  a  rod  bounded  by  an  in¬ 
compressible  shell,  the  compression  of  the  medium  corresponds  to  a 

FTD-TT-63-360/1  +2+4  -2- 


uniaxial  deformed  state.  Thus 


=  -  p,  *  (F  - 1'») '  v„  (1. 1) 

By  virtue  of  (1.1),  the  relation  p  =■•  p(v)  will  be  linear  if  the 
relation  a  =  a(e)  is  linear. 

The  equations  of  the  factors  approximating  the  compression  curve, 
in  a  system  of  units  p  and  V,  are 

/'  ="  —  A*v  "tan  /'  '■>/’>„.  P  -----  —  AJV  ■;  11.,  wbon  P>pm  /j_t  2) 

(.1,  li  finiri.)  '  ’  • 

The  thermal  energy  losses  on  compression  of  the  medium  are"  deter¬ 
mined  by  the  area  of  the  figure  on  the  plane  pV  bounded  by  a  straight 
line  which  connects,  in  the  compression  diagram  p  =  p(V),  the  pointB 
corresponding  to  pressure  at  the  front  and  ahead  of  the  front  and  to 
an  unloading  line.  Therefore,  by  the  first  unloading  law  (dV/dt  «  0), 
the  energy  losses  by  wave  motions  are  as  high  as  possible  with  the 
given  compression  law  p  «*  p(V). 

2.  We  will  examine  the  interaction  between  a  wave  and  a  wall 
when  the  unloading  line  is  parallel  to  axis  £.  Let  us  assume  that 
when  t  =  0  in  the  initial  cross  section  of  the  medium  h  =  0,  the  pres- 
•  sure  jumps  to  p  ,  and  then  drops  according  to  the  given  law 

/»*•=/(/)  (2.1) 


Fig. 


FTD-TT-62-260/1+2+4 


-5- 


The  shock  wave  begins  to  propagate  throughout  the  medium.  The 
flow  of  the  shock  wave  "behind  the  front,  in  Lagrangian  coordinates 
where  h  is  mass  and  t_  is  time,  is  defined  "by  the  motion  equations 


till  , 

7—  '  rr  tai 

01  1  0  4  • 


0, 


Ou  i)V 
oil  ~  oi 


(2.2) 


Here  u  is  the  particle  velocity. 

In  our  case,  unloading  of  the  medium  occurs  behind  the  front 
(region  1  in  Pig.  2)  when  5v/3t  =  0,  from  where  we  obtain  the  solution 
of  Eq.  (2. 2) 

%  =  0.  a  ■-  •  <f,  (/).  I '■  -  <f ,  (/),  '  i>  --  (/)  -i  Mh  (0  (2.3) 


The  functions  cpi(t)  and  ^i(t)  must  be  determined  from  the  initial 
and  boundary  conditions  of  the  problem. 

When  h  =  0,  we  have  p  =  f(t).  Hence  ^i(t)  =  f(t).  The  relations 
at  the  front  of  the  shock  wave  in  the  coordinates  h,  t  have  the  form 

P  —  v o  •"  (Vo  —  V'),  u  —  i/„  -  h  (Vu  —  V)  (2.4) 

where  h  is  the  front  velocity.  We  will  assume  the  pressure  and  veloc¬ 
ity  of  particles  in  front  of  the  wave  to  be  p0  =0,  u0  ®  0.  Let  the 
maximum  pressure  in  the  cross  section  h  =  0  corresponded  to  maximum 
value  pm  in  the  first  part  of  the  approximation.  Then  h  -  An  On 
the  line  of  the  front  h  =  Ait  we  have 


p=Alu,  f  (0  -  /«[),  (()  -  7il(pi  (I)  =  o  (2.5) 

Integrating  the  last  equation,  we  find  cpi(t).  In  this  way  the 
solution  in  region  1  is  obtained 

P  {h,  t)  =  —  (l)  +  /  (i).  u  =>  cpi  (0  (2.6) 


If  the  initial  cross  section  £_  is  given  as 
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MO- 


(2-7) 


then  integrating  (2.5),  we  find  that  in  region  i 


a 


■ii 


...  1 1  M 

{0  •'  7f,  V1  “  afii '  • 


(2.8) 


Let  us  assume  that  a  wall  (or  the  boundary  of  media)  is  situated 
in  cross  section  h  =  h*.  The  relation  p  =  p(V)  in  the  second  medium 
behind  the  wall  is  approximated  by  a  straight  line 


/>  --  -  .r-v  -  :  ir 


■  (2.9) 


When  t*  =  h*/Ai,  the  front  of  the  wave  reaches  the  wall.  Upon 
reflection  form:  region  2  —  reflected  shock  wave,  region  3  —  reflected 
plastic  wave,  region  1*  —  a  passing  wave  behind  the  barrier.'  In  re¬ 
gion  2  loading  of  the  medium  occurB  along  the  unloading  line  (V  «  V(h)) 
up  to  the  limit  of  elasticity  which  is  different  for  different  par¬ 
ticles.  The  solution  in  region  2  relative  to  (2.3)  is 

I'  “  ~  Ht  (0  -i-  '|'s  (/)..  a  --  i;-.,  (/)  (2. 10) 

When  h  =  0,  we  have  p  =  f(t).  Hence  ip 2(t)  «  f(t).  The  front 
velocity  is 

‘ (*•») 

The  velocity  is  determined  by  the  angle  of  Inclination  of  the 
secant  in  the  diagram  p  =  p(V),  drawn  from  a  point  where  p2  corre¬ 
sponds  to  the  pressure  at  the  front  to  a  point  where  pi  is  the  pres¬ 
sure  ahead  of  the  front.  If  the  time  of  the  wave  action  8  in  cross 
section  h  =  0  is  great,  then  the  maximum  pressure  in  the  Incident 
wave  insignificantly  decreases  with  the  increase  of  h.  The  reflected 
wave  moves  over  the  region  with  little  change  of  pressure.  Therefore. 


I 


the  front  velocity  of  a  reflected  wave  (when  p  >  pm)  can  be  assumed 
constant  and,  by  virtue  of  our  approximation,  equals  -Aa.  Then  the 
equation  of  the  line  of  the  front  2-3  (Fig.  2)  is 

h  (2.12) 

On  this  line  in  region  2  the  pressure  in  each  particle  reaches 
the  value  which  was  at  the  front  of  the  Incident  wave  h  =  Ait.  Thus 
by  virtue  of  (2.10) 


t*  (0  =  >f.  M  1/  (/)  -  /  (>.)’■  /  /ill 

1  »  !(/!,  -!-  ,1..)  /•  -  AJ\  I  /.j 


(2.13) 


Integrating  this  equation,  provided  that  q>2(t*)  =*  <Pi(t#),  we 
find  <pa(t)  which  together  with  the  condition  ^/a(t)  =  f(t)  determines 
the  flow  in  region  2. 

If  when  h  »  0  the  pressure  is  defined  by  Eq.  (2.7),  then  from 
(2.13)  we  will  find  that  in  region  2 


*!-  (')  ■ 


I  Mi 


2/1.  L 


<h(l) 


2/i|r  1 1  -l')  3  ( .1  i  /la)  h • 

/iso  1  .Tiao 


/* 

o' 


(2. 14) 


Unloading  of  a  medium  occurs  in  region  3i  the  solution  has  the 
form  of  (2.3)» 

The  condition  at  the  wall:  the  velocity  of  the  particles  of  the 
medium,  adjoining  both  sides  of  the  wall,  equals  the  wall  velocity 
cp3(t).  The  pressure,  acting  from  the  side  of  the  second  medium  when 
condition  (2.9)  is  fulfilled,  is  related  to  the  particle  velocity  by 
the  ratio  p  =  A*u  «  A*q>3(t). 

Thus  the  equation  of  wall  motion  in  region  3  is 

"'<f3  (0  =  -  A*.?,  (!)  +  (/)  -  /l*(p3  (z) 
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where  m  is  the  mass  of  the  wall  per  unit  cross-section  area.  Index  3 
pertains  to  region  3.  We  find  from  the  condition  the  function  ^3(t) 
at  the  front  of  a  reflected  wave.  According  to  (2.4)  and  (2.3) 

Pa  ~Pi  **  —Ai  ("s—  «s).  “  Ha  (0  4-  Hs  (0  +  H'.i  (0  —  /  (0  =  .  . 

—  yla  [tp.,  (/)  —  tp3  (/)]  (2.15) 

Having  eliminated  by  means  of  (2.15)  the  function  ^3,  we  will 
transform  the  equation  of  motion  to 


(A  4  —  A.,1*  “|-  m)  c‘p.t  (/)  =  —  (Ha  -i-  /l*)  q>3  (() 

-!-  H8q>4  (!)  -  \(A,  -I-  At)  t*  -AS  (/)  -|-  /  (/) 


(2.16) 


Integrating  this  equation,  provided  that  q>3(t*)  =  0,  we  will 
find  the  velocity  of  the  wall  q>3(t).  Having  determined  ^s(t)  from 
(2.15),  we  obtain  the  solution  in  region  3.  If  the  wall  Is  rigid 
(m  =  «),  then  in  region  3  the  particle  velocity  is  equal  to  zero 
cp3(t)  »  0.  From  (2.15)  we  will  obtain 

/>(<)-  'M0-  Ms(0-l-/(0-  (2.17) 

If  condition  (2.7)  Is  fulfilled  in  cross  section  h  »  0,  then 
In  region  3  the  pressure  1b  defined  by  the  equation 


Pit) 


ilL±_d»n  (\  1  _  iL'* 

A  i  '’<«  \  "■  ij  Y) 


A,  I*  >.\ 
•la  0  "  t’l 


Hence  we  obtain 


(2.18) 


Aj  jh  ,  l .  r\ 

1  .-I,  Am(l -- j,-0)  when  /-„•=/* 

In  an  Incident  wave  at  that  instant  of  time  according  to  (2.8) 

*  *  hi  \  20- 


Hence  the  coefficient  of  reflection  from  the  wall  is 
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(2.19) 


If  when  h  a  0  the  pressure  is  given  in  the  form  of  (2.7)  and  the 
mass  m  is  finite,  then  introducing  a  new  variable  t  =>  t  -  t*,  we  will 
obtain  the  differential  equation  of  wall  motion  (2.16)  in  region  5  in 
the  form 


("<  A.r)  <i>:,  (t)  (/ls  4-  A*)  <[>„  (r)  •• 

-U  N;  JL il  iilil  1,1  (l .. 

.-I,  L  •  0  20  1  .1.0 


1,1 


(2.20) 


For  large  values  of  t*/0  the  pressure  drop  at  the  front  of  the 
incident  wave  in  the  section  0  <  h  <  h*  is  great  and  the  front  veloc¬ 
ity  of  a  reflected  plastic  wave  can  be  sharply  distinguished  from  -As- 
In  this  case  the  flow  in  regions  2  and  3  must  be  obtained  by  the  simul¬ 
taneous  determination  of  the  boundary  between  them. 

If  the  relation  p(V)  on  compression  is  approximated  by  one 
straight  line  having  an  inclination  Ai  and  unloading  occurs  along  the 
vertical,  then,  assuming  in  Eq.  (2.20)  A*  =»  kz  =*  Ai,  we  will  obtain 
its  solution  as 


,  .  ^i.i  u .  r  m  \  ,  2  in  /  I  ..  \  . 

W  “  X  U1  --  af --  TiTu)  ( 1  ~ **)*'•  3  -io  la  - a‘  j  ••• 
TTlY lA*  a)  ~v  •'  t)  I  +  V  -  Wl 1,1  jr=a7?|| 


Here 
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The  front  velocity  of  a  reflected  plastic  wave  is  equal  to  -Ai. 

Depending  on  the  parameters  of  the  incident  wave,  on  the  charac¬ 
teristics  of  the  media  on  both  sides  of  the  wall,  on  the  mass,  and  on 
the  removal  of  a  wall  from  the  first  cross  section,  various  configura¬ 
tions  of  a  reflected-wave  system  are  possible  (in  particular,  unload¬ 
ing  regions  can  be  replaced  by  loading  regions,  etc.  ). 

3.  We  will  examine  the  interaction  between  a  wave  and  a  wall  in 
the  case  where  the  unloading  line  is  not  parallel  to  axis  Let  us 
assume  in  cross  section  h  =  0  that  the  pressure  changes  in  conformity 
with  (2.7).  We  will  approximate  the  curve  of  compression  of  the  med¬ 
ium  by  a  straight  line  (1.2),  and  the  unloading  line  curve  by  a 
straight  line 

P  —  —  -'s8W  -I  -  II *  _  (3.1) 

where  B2  depends  on  the  value  of  maximum  stress  attained  on  compres¬ 
sion;  pm  corresponds  to  maximum  pressure  in  the  first  section  of  the 
approximation.  The  solution  of  Eqs.  (2.2)  defining  the  flow  behind 
the  wavefront  formed  in  the  medium,  for  linear  approximation  of  the 
relation  p  =  p(V),  as  was  shown  in  [l]  is 

/'  i'\  ('*  -  'UO  "!*  (/*  +  •»„/),  .\,u  »  I'\  (h  -  ,u).  -  (h  +  ,1,/)  (3.  2 ) 

The  functions  Fi  and  F2  are  defined  by  the  Initial  and  boundary 
conditions.  The  flow  behind  the  wavefront  (region  1  in  Fig.  3) >  ob¬ 
tained  from  (3-2),  according  to  [l],  is  defined  by  the  equations 


I,  Ar-'-.lrh  »\ 

P  Pm\l  'ii.l,.!,*  “if  "0  i 

*  p,n  i'  li  '  to  iMi.  i,  ?) 


(?-3) 


When  t  =  hx/Ax  *  t  ,  the  front  reaches  the  wall.  In  this  case 
regions  2  and  lx  are  formed.  As  In  part  2,  we  will  assume  the  equatic 


o_ 


of  the  reflected  wavefront  in  the  form  of  (2.12)  and  the  relation 
P  =  p(V)  in  the  second  medium  behind  the  wall  in  the  form  of  (2.9). 


The  region  corresponding  to  region  2  in  Fig.  2  in  the  case  con¬ 
sidered  does  not  appear,  since  the  velocity  of  the  reflected  wave- 
front  coincides  with  the  unloading  velocity.  The  condition  at  the 
boundary  of  1-2  (Fig.  5) 

Pz  +  Atu.i  =  +  /11«1 

yields  that  the  function 


transfers  from  region  1  to  region  2. 

Thus  the  differential  equation  of  wall  motion  is 


(3.*) 


(3-5) 


or 


where 


««F  (0  =  Pi  —  PC'  21' i  (/<y  —  a ../)  —  (A,  *!■  A  ,>)  <p  (0 


<i-  -i-  <>  !-  m  -i-  u  —  o 


/ >  .. 


■  i,  -  ij 
P’»  Ml.  I, 


(A i  -I-  ✓l-Wi’* 
2At-\l 


(3.6) 


2.1,.  I,  0m  ’ 

Integrating  (5.6),  provided  that  condition  <p(tx)  =0,  we  will 
find  the  wall  velocity 


•i- 


it  f  ,  ,  /,x  .  /  i  i)  i,'\  (.  ,  /,  h*\\\ 

w  ■“  it  {” '  -■  nr  (c  cxp  [-  c  ('  -  ~j\)} 


(3-7) 


We  will  determine  the  second  function  in  region  2  from  the  con¬ 
dition 
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U)  -  .l,.r.  (/) 


Hence 


F,  (If  -  AJ)  -  F,  (If  . 


(h  +  A tt) 


■  U  -!-  .  ia 
2,1! 


■  if + /ta(\ 
At  1 


The  solution  in  region  2  is 


i‘  - 


A, 


A- A-.  |,  -li-r.lj,!  ,  III— If- 1- 

AT-l,A.i'r~tm‘n  "  "rgr'  -  A#i — ,ir~) 


•  I  i  -1-  At 


.•la 


(3.8) 


We  obtain  the  solution  in  region  lx  from  the  condition  at  the 
front  of  the  passing  wave  px  -  A^i  =>  2FX  =  0  and  from  the  condition 
that  when  h  =  hx  the  particle  velocity  equals  the  wall  velocity 

/'  -  A*u  . l-\«  (h  -  A-t)  --  ,i>  l 


Pig.  4. 


The  wall  velocity,  defined  by  equation  (3*7  )>  first  increases 
and  then  drops  off.  Differentiating  (3.7)  and  equating  the  derivative 
to  zero,  we  will  find  that  the  maximum  value  of  velocity  is  reached 
at  the  instant  of  time  t,  determined  from  the  condition 


.1'  ■  I  •.!.!  .0  -r  —  .l.)>ix)  V * 

.  Ij  ■ .  I  it  i  1, .  I  >in  j 


(3-9) 


If  m  =»  »,  then  in  region  2  when  h  «=  h? 


«  «  0,  F ,  (A*  -  .1,0  -  F,  (A*  *1*  >1.0  -  0 


Hence 


/•'=  (*  -i  .1.0  -  (2A*  -  A  -  ^,i) .  />m  X 

w‘<*  •  •  <*4  *- 


(5.10) 


If  m  =  0,  then  from  (3*6)  we  will  find  that  when  h  =  hx 


<p(0 


(i>- -  .1,0 
F4-  ,!, 


(- 1.  - .  Ia)  /',» 
(.1-  .  iT>TT7  1 


.ii  f-  .i. 


(l<  - 


(3.11) 


4.  We  will  analyze  the  regularity  of  wall  movement  in  an  elemen¬ 
tary  case  where  the  relation  o(s)  on  compression  is  approximated  by 
one  straight  line. 

We  will  examine  the  calculation  results  of  the  wall  velocity, 
which  are  presented  in  Fig.  4  by  two  graphs  corresponding  to  the  dif¬ 
ferent  models  of  the  medium. 

In  both  cases  the  wall  mass  m,  its  distance  from  the  initial 
cross  section  hx,  and  the  pressure  change  in  cross  section  h  =  0 
(Eq.  (2.7))  are  assumed  identical. 

The  graphs  correspond  to  the  same  law  of  compression  of  a  medium 
(one  linear  factor),  however,  the  unloading  line  is  different:  curve 
1  is  the  unloading  line  coinciding  with  the  load  line;  curve  2  is  the 
unloading  line  parallel  to  axis 

In  the  second  case  under  these  conditions,  the  energy  loss  on 
movement  of  the  incident  and  reflected  waves,  as  noted  above,  are  max¬ 
imum,  and  in  the  first  case  are  minimum. 
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Fig.  5a. 


Fig.  5"b. 


The  difference  in  the  character  of  the  graphs  is  due  to  the  dif¬ 
ference  in  energy  loss.  The  greater  the  loss,  the  smaller  the  value 
of  maximum  wall  velocity. 

Figures  5a  and  5b  show  the  pressure  change  at  the  wall  -during  its 
movement  for  these  two  models  of  the  media:  Fig.  5a  is  the  medium 
governed  by  Hooke's  law;  Fig.  5b,  the  unloading  line  is  vertical.  In 
Fig.  5a  and  5b  the  pressure  in  the  wave  In  cross  section  h  =»  hx  In  the 
absence  of  a  wall  Is  shown  by  a  dashed  line.  It  is  obvious  from  the 
figures  that  the  presence  of  a  wall  appreciably  distorts  the  wave  only 
In  a  time  interval  small  as  compared  with  0. 

Thus  in  an  elasto-plastic  medium  with  vertical  unloading.  Just 
as  in  a  linearly  elastic  medium  governed  by  Hooke's  law,  even  a  wall 
with  a  relatively  large  nass  is  involved  in  motion  together  with  the 
medium  during  a  time  interval  small  as  compared  to  the  time  of  the 
wave  action. 

This  conclusion  also  pertains  to  media  with  an  Inclined  unload¬ 
ing  line. 
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